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ABSTRACT 

Preliminary  experiments  in  the  numerical  prediction  of  large  scale 
sea-surface  temperature  anomalies  are  made  using  a  10-level  primitive 
equation  model  with  300-km  horizontal  resolution  covering  a  rectangular 
basin  in  the  North  Pacific.   The  model  is  first  integrated  over  an  11- 
year  period  to  statistical  equilibrium  using  time  dependent  wind  and 
thermal  forcing.   The  monthly  normal  climatology,  generated  by  the 
model  in  this  long  term  integration,  is  then  used  along  with  observed 
sea-surface  temperature  anomalies  to  define  the  initial  state. 

The  results  show  that  inclusion  of  salinity  in  the  model  may  pro- 
duce more  accurate  predictions  of  high  latitude  cold  anomalies.   The 
most  accurate  of  the  predictions  removed  all  of  the  anomalous  diffusion 
terms  in  the  governing  equations.   Data  for  the  observed  temperature 
anomalies  were  not  available  at  depth,  thus  initialization  of  the 
anomalies  below  the  surface  was  superficial.   Consequently,  the  effects 
of  vertical  advection  were  not  fully  realized,  indicating  a  need  for 
future  experimentation  with  a  more  accurately  defined  initial  profile 
of  the  anomalies  in  the  initial  conditions. 
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I.   INTRODUCTION 

For  today's  meteorologist  to  make  accurate,  long  range  forecasts  the 
understanding  of  the  dynamics  of  atmospheric  behavior  is  imperative.   In 
the  recent  meteorological  past,  the  importance  of  air-sea  interactions 
has  been  realized  and  has  given  the  forecaster  another  potential  tool  with 
which  to  improve  his  predictions.   A  better  understanding  of  the  inter- 
actions between  the  oceans  and  the  atmosphere  will  also  assist  in  ocean 
acoustic  forecasting  and  other  Naval  operations. 

Like  the  atmosphere,  the  oceans  experience  seasonal,  monthly,  and 
daily  variations  in  temperature  and  circulation  patterns.   Further  investi- 
gation of  these  variations  reveals  the  existence  of  anomalous  behavior 
within  the  mean  temperature  and  circulation  patterns.   It  is  believed 
that  these  variations  influence  directly  the  atmospheric  weather  patterns 
through  the  exchange  of  moisture,  momentum,  and  heat  across  the  air-sea 
interface.   This  thesis  is  concerned  with  the  prediction  of  these 
anomalous  oceanic  temperatures  and  currents. 

During  the  past  two  decades  it  has  become  apparent  through  the  studies 
of  J.  Namias  and  the  late  J.  Bjerknes  that  the  atmospheric  response  to 
sea-surface  temperature  anomalies  may  be  worldwide.   Bjerknes  (1966,  1969, 
1972)  has  shown  that  large-scale,  sea-surface  temperature  anomalies  exist 
in  the  equatorial  region  and  that  these  anomalies  produce  a  response  in 
the  atmosphere  which  is  not  limited  to  that  region  alone.   Likewise, 
Namias  (1959,  1969,  1972b)  discusses  the  possibility  that  such  a  response 
is  also  present  in  the  middle  latitudes.   Furthermore,  Namias  (1969) 
observed  that  during  several  years  when  intense  sea-surface  temperature 


anomalies  were  present  in  the  North  Pacific  Ocean,  there  were  distinct 
weather  patterns  over  the  United  States.   Because  of  the  above  pioneer- 
ing findings  of  Namias  and  Bjerknes,  the  North  Pacific  Experiment 
(NORPAX)  was  established  in  which  a  cooperative  scientific  effort  is 
being  made  to  understand  the  causes  of  these  anomalies  and  their  rela- 
tionship to  the  atmosphere. 

Numerical  integrations  of  coupled  ocean-atmosphere  models  have  been 
reported  by  Manabe  (1969) ,  Bryan  (1969) ,  and  Wetherald  and  Manabe  (1972) 
By  comparing  the  results  of  their  coupled  model  with  those  of  the  ocean 
and  atmosphere  models  integrated  separately,  these  authors  showed  the 
significant  role  ocean  circulation  plays  in  determining  the  state  of 
the  atmosphere. 

In  turning  our  attention  to  the  oceans,  the  first  attempts  at  large 
scale  ocean  predictions  was  based  upon  the  use  of  a  prediction  equation 
for  the  sea-surface  temperature  (Namias,  1965).   In  Namias1  method, 
temperature  T  is  a  function  of  x,y,z,  and  time  t.   The  local  rate  of 
change  of  temperature  can  be  represented  by 


dl               8T     8T     3T 
■kt-  -  -U  -k v-r w^—  +q  (1) 

dt      dx     dy     dz 


where  q  is  the  rate  of  change  of  temperature  from  non-advective  pro- 
cesses.  Now  if  we  let  (  )  denote  the  long  term  normal  state  and  (')  the 
anomaly,  then  the  monthly  mean  can  be  expressed  by 

T=T+T'  ,u=u+u'  ,  v  =  v  +  v1  ,w=w+w' 

and 

q  =  q  +  q«  . 
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Then  Eq.  (1)  takes  the  form: 


^-•)=-<G+u-,-ii|^l-(v+v,i<|p 


-(w+w«)   ^  (^T,)  +  i  +  q'  .  (2) 


Assuming  the  long  term  normal  temperature  satisfies  the  equation 


3T     -  3T   -  3t   -  3t 

■5—  =  -  u  -s V  "s w-r—  +q  (3) 

dt       3x     3y      dz 


subtracting  (3)  from  (2)  gives  the  equation  for  the  temperature  anomaly: 


3T'  .,-  .,  3T'  -  ,,  3T'  -  ,.  3t.' 
■5-  =  -  u+u')  •5-+(v+v')  -T-  +  (W+W1)  -5-] 
3t  3x         dy  3z 

r  1  3T  ,   ,  3T  j   ,  3t.     .  ... 

-u'T-+v,-5-+wl3-]+q'.  (4) 

3x      3y      3z 


Namias '  first  predictions  (1965)  retained  only  the  advection  of  the 

3t      3t 
normal  temperature  by  the  anomalous  currents,  (u'  v—  +  v1  ■r— )  ,  in  Eq. 

(4).   The  anomalous  currents  (u* ,  v1)  were  obtained  from  the  anomalous 

atmospheric  wind  stress  observed  during  the  period  of  the  "forecast" 
and  the  long-term  normal  temperatures  (T)  from  observations  collected 
by  the  U.  S.  Weather  Bureau.   Since  this  prediction  used  data  actually 
observed  during  the  forecast  period,  it  is  not  a  prediction  in  the  true 
sense  and  is  referred  to  as  a  "semiprediction" .   Although  this  simple 
method  produced  results  which  correlated  positively  with  the  observed 
temperature  changes,  Arthur  (1966)  demonstrated  that  additional  terms 
should  be  included  in  the  prediction  equation.   As  a  result,  Jacob  (1967) 
performed  a  semiprediction  neglecting  only  the  vertical  advection  terms 
and  the  anomalous  heating.   A  quantitative  comparison  of  his  results 
indicated  an  improvement  in  the  prediction. 
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More  recently,  Namias  (1972a)  has  shown  that  the  sea-surface  tempera- 
ture patterns  in  the  North  Pacific  can  be  predicted  with  some  success 
using  normal  surface  currents  to  advect  the  anomalies  (neglecting  all  but 

—  3T '    —  r)T ' 
-u  s—  »  -v  -s—  on  the  right  side  of  Eq.  (4)).   The  results  of  his  study 
dx      dy 

(Figs.  1,  2,  3,  and  4),  which  will  be  discussed  in  detail  later,  show 
the  importance  of  horizontal  advection.   These  results  plus  the  lack  of 
temperature  data  at  depth  clearly  indicates  the  need  for  a  completely 
dynamical  model.   Therefore,  it  is  the  objective  of  this  thesis  to  predict 
the  evolution  of  the  anomalous  temperatures  using  a  numerical  model  which 
contains  all  the  advection  terms  in  (4)  and  a  parameterization  of  mixed 
layer  dynamics  as  well. 

To  make  a  comparative  analysis,  the  period  covered  by  this  prediction 
will  be  the  same  as  that  used  by  Namias  (1972a) .   Sea-surface  temperature 
anomaly  data  were  also  obtained  from  Namias'  20-year  normal  observations. 
It  should  be  noted,  however,  that  this  normal  data  is  from  the  years  1947 
to  1966  and  is  not  necessarily  representative  of  the  prediction  period. 
Furthermore,  due  to  the  aforementioned  lack  of  data  at  depth,  we  assume 
an  anomalous  temperature  profile  which  decreases  linearly  from  its  value 
at  the  surface  to  zero  at  150  meters.   NORPAX  scientists  have  shown  that 
anomalies  exist  at  depths  greater  than  150  meters;  however,  this  depth 
was  chosen  to  correspond  more  closely  to  Namias1  study  which  assumed  the 
temperature  anomaly  to  be  confined  strictly  to  the  surface.   Further 
studies  could  be  conducted  varying  the  vertical  profile  of  the  tempera- 
ture anomaly  to  determine  the  importance  of  the  initial  condition.   The 
initial  anomalous  current  is  geostrophic  and  its  vertical  average  is 
zero.   Unlike  Namias'  method  which  employed  the  use  of  observed 


12 


climatological  data  during  the  entire  forecast  period,  this  method  uses 
model-generated  climatology  only  to  define  the  initial  conditions.   Con- 
sequently, following  initialization  the  anomalous  and  mean  fields  are 
completely  free  to  dynamically  interact. 
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II.   THE  NUMERICAL  MODEL 

A.   BACKGROUND 

The  basic  model  used  in  the  following  predictions  is  that  of  Haney 
(1974)  as  modified  by  Davies  (1975) .   "The  model  is  based  upon  the  hydro- 
static and  Boussinesq  approximations.   This  means  that  the  ocean  is 
assumed  incompressible,  and  the  density  is  replaced  by  a  constant  every- 
where except  in  the  hydrostatic  equation  where  it  is  multiplied  by  the 
acceleration  of  gravity.   A  major  simplification  is  the  neglect  of 
salinity  and  the  assumption  that  the  density  is  a  linear  function  of 
temperature  alone." 

The  equations  which  govern  Haney ' s  model  are  given  by 

du     -1    3p. _  u         ,    „  „2      92u 

—  = -ttF  +fv  + v  sin  0  +  A  V  u  +  K ,             (5) 

dt    p  acosd)  dk  acos©            m         „  2 

o  dz 

2 

dv    -1  3p  r      u      .         „2      9  v 

—  =  -^  -fu r  u  sin  4>  +  A  V  v  +  k ,  (6) 

dt   p  a  del)       acosi)  m         »  2 

o  dz 

|f  -  -pg  ,  (7) 

3w     1   r9u    3  ,      , .  . 

■5-  +  r  TT  +  ^T  v  cos  4>)  1  =  0  f  (8 

dz   acostp  dX        6$ 

2 

—  =  A   V  T  +  K +  6  (T)  ,  (9) 

dt     H  «  2     c 

dz 

p   =  p   [1  -  a(T  -  T  )]  ,  (10) 

o  o 

where  the  representation  of  the  symbols  is  as  follows: 

t   time 

X    longitude 

cf>   latitude 


Haney,  R.  L. ,  "A  Numerical  Study  of  the  Response  of  an  Idealized 
Ocean  to  Large  Scale  Heat  and  Momentum  Flux,"  J.  of  Physical  Oceanog., 
V.  4,  No.  2,  p.  147,  1974. 
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z  height 

ft  angular  speed  of  the  earth's  rotation 

f  coriolis  parameter 

a  radius  of  the  earth 

g  acceleration  of  gravity 

u  eastward  velocity  component 

v  northward  velocity  component 

w  vertical  velocity  component 

T  temperature 

T  constant  reference  temperature 

o 

p  density 

p  density  of  water  at  reference  temperature 

a  coefficient  of  thermal  expansion 


B.   BOUNDARY  CONDITIONS 

After  Haney  (1971)  and  Davies  (1975) ,  the  determination  of  the  upper 
boundary  condition  in  the  thermal  energy  equation  is  computed  by  calcu- 
lating separately  the  incoming  solar  radiation  (Q  ) ,  the  net  upward  long- 
wave radiation  (Q  ) ,  the  latent  heat  exchange  (Q  ) ,  and  the  sensible  heat 
B  E 

exchange  (Q  ) .   Together,  these  fluxes  define  the  net  downward  heat  flux 
H 

as  represented  by 


2  =  QI    ~  (QB  +  QE  +  V  '  (11) 

The  back  radiation  Q   is  calculated  from 

*B 

Q„  (T  )  =  .985(.39  -  .05e  1/2) (1  -  .6C2)  0   T  4  (12) 

Bo  a  o 

where  T   is  the  sea  surface  temperature,   e   the  atmospheric  vapor 
pressure  at  10  meters,   C   the  fraction  of  cloud  cover,  and  a  the 
Stefan-Boltzmann  constant.   The  exchange  of  latent  heat  and  sensible  heat 
across  the  ocean-atmosphere  interface  were  computed  from 


Qv   =    .622  p   C^  V  L  [e  (T  )  -  e  ]/P  (13) 

E        a  D  '  ■     so     a   a 
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and 


Q   =  p   C    V   C   (T   -  T  )  , 
*H    Ka  D  '  '   p   o    a 

respectively,  where  p   was  the  density  of  the  air,   C   the  drag  co- 

efficient,   |v|   the  wind  speed  at  anemometer  level,   L  the  latent  heat 

of  vaporization,   P   the  atmospheric  pressure,   e  (T  )   the  saturation 

a  so 

vapor  pressure  at  the  sea-surface  temperature,   C   the  specific  heat 

XT 

of  air  at  constant  pressure,  and  T   the  atmospheric  temperature  at 

Si 

10  meters.   Incoming  solar  radiation  (Q  )  was  calculated  with 

QT  =  S  (.74  -  .6C)  (15) 

I    o 

where  S   was  the  solar  insolation  incident  on  a  horizontal  surface  at 
o 

the  top  of  the  atmosphere. 

The  quantities   S   ,  C  ,  e   ,  T   and  |v|  are  prescribed  functions 

of  latitude  (v  depends  on  longitude — see  below)  and  time  taken  from 

atmospheric  climatological  data.   The  fluxes  Q   ,  Q   ,  and  Q   are 

BE         H 

additionally  dependent  on  the  computed  sea-surface  temperature,  T 
Therefore,  an  anomalous  sea-surface  temperature  produces  a  corresponding 
anomalous  heat  flux  which  has  the  tendency  to  remove  the  temperature 
anomaly  by  virtue  of  the  opposing  sign  in  Eq.  (11) . 

"In  addition  to  the  thermal  forcing,  the  ocean  circulation  was  also 
forced  by  surface  winds  given  on  a  time-dependent  basis.   Assuming  an 
average  geostrophic  inflow  angle  of  10°  and  a  frictional  reduction  of 
velocity  by  10%,  one  could  derive  the  surface  stress  from  the  surface 
geostrophic  winds.   If  u  and  v  are  the  eastward  and  northward  surface 
winds,  respectively,  after  inflow  and  frictional  influences  are  considered, 

then 

ii      2     2  1/2 
V  =  (u  +  v  )A/ 
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and  the  wind  stress  components  are 

T   =  p  C   I V I  u 

x    a   D  '  ' 

I  l     ,.2 
T   =  p   C   V  v  ." 

y    a  D  '  ' 

Therefore,  the  wind  stresses  are  completely  specified  since  u(y,t) 
and  v(x,t)  are  completely  independent  of  the  ocean  temperatures  and 
currents  in  this  experiment.   Since  the  quantities  which  define  the 
atmospheric  forcing  are  derived  from  climatological  data,  this  thesis 
does  not  concern  itself  with  the  response  of  the  ocean  to  anomalous 
atmospheric  forcing. 

Lastly,  the  ocean  bottom  is  assumed  to  be  flat  which  requires  the 
vertical  velocity  to  vanish  there.   Thus  integrating  the  continuity 
equation  up  to  the  surface  and  applying  a  free  surface  boundary  condi- 
tion, one  could  obtain  a  prediction  for  the  height  of  the  surface. 
However,  the  present  method  assumes  that  the  ocean  surface  is  not  a  free 
surface  but  rather  a  balanced  surface  at  which  w  =  0.   Therefore,  since 
w  =  0  at  the  top  and  bottom,  the  vertical  average  of  (8)  shows  that  the 
vertical  average  current  is  non-divergent  and  consequently  can  be  de- 
scribed by  a  streamfunction.   The  streamfunction  is  predicted  from  the 
vorticity  equation  which  is  derived  from  the  vertical  average  of  the 
equations  of  motion. 

C.   RESOLUTION  AND  DATA 

To  improve  the  resolution,  the  6-level  model  employed  by  Haney  (1974) 
was  expanded  to  include  10  levels  (Davies,  1975).   The  vertically  expanded 


2    .  .    . 

Davies,  R.  W. ,  A  Numerical  Parameterization  of  Wind  Mixing  in  a 

Tine-Dependent  Baroclinic  Oceanic  General  Circulation  Model,  M.  S.  Thesis, 
Naval  Postgraduate  School,  Monterey,  California,  1975,  pp  56. 
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model  now  has  fixed  levels  at  10,  30,  60,  100,  150,  225,  350,  700,  1500, 
and  3000  meters. 

As  mentioned  previously,  to  compare  the  results  with  those  of  Namias 
(1972a) ,  sea-surface  temperature  anomaly  data  covering  the  same  period 
were  obtained  from  J.  Namias  of  Scripps  Institution  of  Oceanography. 
The  surface  temperature  data  thus  obtained  covered  only  the  area  20°N  to 
60°N,  130°E  to  110°W  in  the  North  Pacific  and  the  grid  spacing  intervals 
for  the  latitude  and  longitude  were  5°  and  5°  respectively.   Therefore, 
Namias1  data  were  expanded  to  the  equator  and  to  65°N  by  extrapolating 
the  anomalous  data  to  zero  at  the  most  distant  boundaries.   Additionally, 
the  data  over  the  entire  domain  were  interpolated  to  fit  the  2.0°  lati- 
tude and  the  2.8°  longitude  grid  spacing  of  the  present  model. 

As  a  result  of  the  expansion,  the  predictions  outside  the  original 
data  field  should  be  ignored.   Consequently,  for  the  convenience  of  the 
reader,  an  overlay  (Fig.  5)  with  the  locations  of  ocean  vessel  stations 
November  (N) ,  Papa  (P) ,  and  Victor  (V)  as  well  as  the  site  of  the  NORPAX 
"Pole"  experiment,  has  been  superimposed  on  the  normal  fields  and  each 
related  prediction.   The  superficial  land  mass  boundaries  of  the  NORPAX 
grid  are  represented  by  the  shaded  area  and  the  original  data  field  by 
the  area  between  the  dashed  lines. 

D.   INITIAL  CONDITIONS 

To  completely  define  the  initial  state  of  the  ocean  for  dynamical  pre- 
dictions, a  considerable  amount  of  model  climatology  was  used  along  with 
the  available  observed  data.   The  model  climatology  was  generated  by  in- 
tegrating the  model  over  an  11-year  period  starting  from  a  stratified 
initial  state  and  using  atmospheric  data  which  contained  the  observed 
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normal  seasonal  data.   This  model  climatology  was  used  for  initializa- 
tion and  verification  as  described  below.   To  provide  continuity  in 
the  atmosphere  wind  and  thermal  forcing,  the  annual  mean  and  amplitude 
and  phase  of  the  first  two  harmonics  were  calculated  from  atmospheric 
climatological  data  as  defined  by  Davies  (1975)  .   The  atmospheric  data 
were  averaged  along  each  grid  latitudinal  line,  except  for  v,  the  north- 
ward component  of  the  surface  geostrophic  wind,  which  was  averaged 
along  each  meridian. 

The  resulting  initial  normal  currents  and  the  30,  60,  90  day  normal 
currents  at  level  2  (30  m)  are  shown  in  Figs.  6,  8,  10,  and  12,  where 
"initial"  is  defined  as  October  15th  of  year  11.   As  seen  in  Figs.  6  and 
12,  there  is  an  almost  negligible  change  in  the  normal  currents  during 
that  period.   The  corresponding  initial,  30,  60,  and  90  day  normal  tem- 
perature fields  are  shown  in  Figs.  7,  9,  11,  and  13.   Although  there  is 
little  change  in  the  temperature  gradient  in  the  middle  latitudes,  the 
southward  movement  of  the  main  baroclinic  zone  indicates  a  cooling  trend 
which  was  expected  during  the  winter  months. 

To  begin  the  predictions,  the  total  field  must  be  defined  in  the 
initial  conditions.   Thus,  the  normal  part  of  the  initial  conditions  is 
obtained  from  the  model-generated  climatology  while  the  anomalous  part 
is  obtained  from  observations.   The  observed  sea-surface  temperature 
anomalies  are  extrapolated  linearly  from  observed  values  at  the  surface 
to  zero  at  level  5  (150  meters) .   This  rather  arbitrary  procedure  of 
defining  the  anomalous  temperatures  below  the  surface  is  done  because 
of  the  lack  of  sub-surface  temperature  data. 
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The  total  current  (normal  plus  anomaly)  can  be  represented  by 

/s      * 

w  -  w  +  w  (17) 

i  r° 

where   (  )   is  the  vertical  average  current,  [—  /  (  )dz] , (  )  is  the 

H  J 

— H 
vertical  shear  current,  [(  )-C)\,    and  H  is 

the  total  ocean  depth.   Expanding  (17)  into  its  normal  and  anomalous 

parts  one  obtains 

_        *   A   *    * 

w  =  w  +  w'  =  w  +  \v,+  w  +  w'  <18) 

*     * 

where  W  and  W  are  the  normal  part  of  the  vertical  average  current 
and  the  normal  part  of  the  vertical  shear  current  respectively.   The 

anomalous  part  of  the  vertical  average  current  (W')  is  assumed  to  be 

* 
zero  and  the  anomalous  part  of  the  vertical  shear  current  (\V')  is  the 

geostrophic  current  calculated  from  the  above  temperature  anomaly  using 

the  thermal  wind  relation.   Since 

W  =  -^r  Ik  x  Vp  , 

g     pQf 

taking  the  partial  derivative  with  respect  to  z  produces 

3\V 

Ik  x  V(-pq)  . 


dz  p  f 

o 


Substituting  (10)  for   p 


p  =-p  a  T  +  constant  , 
o 


where  the  gradient  of   p   becomes 


Vp  =  -p  a  g  V  T  . 
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Therefore, 
3\V 


=  -—    Jk  x  VT  (thermal  wind  relation) . 

dz  f 

It  is  now  assumed  that  the  anomalous  part  of  the  vertical  shear  current 
satisfies  this  relationship, 

$■'  =  sa  ,k  x  VT.  . 

dz     f 

Then  integrating  from  -H  to  some  depth   z 

/z 
T'd£  (19) 

-H 

*  * 

where  W*  (-h)   is  unknown.   If  the  vertical  average  of  W1   is  taken, 

then  this  quantity  can  be  set  equal  to  zero  by  definition.   Thus,  from 

(19) 

*  o  o  z 

W«    =  jjj-   J     \*'  (-H)dz   +^-    IkxV^y     (    f     T'dgdz    E    0    , 
-H  -H         -H 

or 

o  z 

\V'  (-H)    +22-    Ik  x  V  i  )     if  T'dOdz   =   0    .  (20) 

* 
Then  (20)  can  be  solved  for  \y'(-H)   and  substituting  into  (19)  to  give 

the  result 


z  o     z 

-H 


W'(z)  =22  Ik  x  V  [J        T'dz  ~  Yi      f     {      f    T'dC)dz]  . 
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In  Eq.  (21)   f   is  not  allowed  to  become  smaller  than  its  value  at 
13.25°N.   This  geos trophic  approximation  to  the  vertical  shear  current 
anomaly  is  then  used  as  the  initial  condition  in  the  present  model. 

Having  defined  the  initial  temperature  and  currents  by  combining 
the  model  generated  normal  fields  with  the  observed  anomaly  fields,  the 
model  is  run  for  90  days  of  simulated  time.   During  this  time,  the  total 
temperatures  and  currents  interact  and  evolve  together.   At  any  given 
time  later  (e.g.,  60  days)  the  predicted  anomaly  is  defined  to  be  the 
total  field  predicted  for  that  time  minus  the  model  generated  normal 
for  that  same  time. 

E.   THE  PREDICTION  EQUATION 

Since  several  cases  were  studied  using  variations  of  the  thermodyna- 
mic prediction  Eq.  (9)  which  governs  the  total  temperature  field,  a 
detailed  analysis  of  that  equation  follows. 

Neglecting  the  convective  adjustment  (6  (T) )  in  the  analysis,  (9) 
can  be  written 

dT   „  „  2       82T  ,_, 

_  =  A  Vh   T  +  K  —  .  (22) 

dz 


Expanding  (22) 

3tj    9t    9t    9t      82t  92t,  92t  ,,_, 

dt      dx     dy     dz       _  2  ,.  2  ,.2 

dx  dy  dz 


22 


and  substituting  further  for  the  normal  and  anomalous  parts,  (23)  takes 
the  form 


9t   3t'    r-  9t   -  3t*   ,  9t    „  9t'  -  3t   -  3T' 
9t   9t        9x     9x      9x      dx     dy     9y 

,  9t    ,  3T'  -  3t   -  3t«   ,  9t    ,  3t'. 

+  V'  7T-+V'  -r—  +  W  —  +W-r—  +  w'  -r—  +  W*  -r—  J 
dy       dy      dz      dz       dz       dz 

,r92T   92T'    32T   92T'.     r92T   92TV  ,_„, 

+  A[ + + +  j  ]  +  K[ 2    2  J  '  (   J 

9x    9x     9y    9y        9z    9z 


Setting  T ' ,  u',  v',  and  w*   equal  to  zero,  then  the  normal  field  can  be 
represented  by 

9t     -  9t   -  3T   -  3t.   ,r32T   92T.   ^  92t  ,__, 

dz        dx      dy      dz         2      2         2 

9x     dy        dz 

Subtracting  (25)  from  (24) ,  the  equation  governing  the  anomalous  tempera- 
tures becomes, 

9t'     -  9t'    ,9t    ,9t'   -  9t'    ,  3t    ,9t'   -  9t' 
9t        9x      9x     9x      dy      dy     dy      dz 

2      2         2 

,3T  t   ,3T»]  +  A  [-2-i   +  -2-4  ]  +  k  -2-4   .  (26) 

+  wSr—  +  w1-^—        -  2     „  2        ^2 
9z     9z       9x     9y        dz 

The  following  investigation  will  reveal  that  initially  \v,#VT"  =  0 
by  analysis  of  (21).   If  we  let  "D"  represent  the  depth  at  which  the  tem- 
perature anomaly  becomes  zero  as  defined  earlier,  then  the  temperature 
anomaly  (T1)  at  some  level  z  <  0  is 


T' (x,y,0)  -  T(x,y,-D)   =   T' (x,y,0)  -  T' 
D  -z 
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where  T' (x,y,0)  is  the  sea-surface  temperature  anomaly  and  T' (x,y,-D) 
is  the  temperature  anomaly  at  D  =  150  meters.  Since  T1 (x,y,-D)  =  0  by 
definition,  then 

T'  =  T'  (x,y,0)  (1  +  §)  .  (27) 

Substituting  (27)  into  term  (a)  of  (21)  gives 

z 

I     T'dz  =  F(z)  T'(x,y,0),  (28) 

-H 

where 

z 
F(z)  =  J      (1  +  |)d? 

-H 

and  (£)  is  a  dummy  variable  representing  the  vertical  coordinate.  Like- 
wise, (28)  can  be  substituted  into  term  (b)  of  (21)  to  produce 

o    z  o 


h  f  (  /  Trd^)dz  =  h  /  F(z)  T'  (x'y'0)dz 


-H    -H  -H 


=  F  T'  (x,y,0)  .  (29) 


Therefore,  (21)  becomes 


*    ccq      „  A 

W'  =  -j-    Ik  x  VH[F(z)  T'  (x,y,0)  -FT'  (x,y,0)] 


or  by   factoring      T" (x,y,0) 


*  ota  « 

W    =  -f-    Ik   x  VfT"  (x,y,0)     (F(z)    -   F)]     . 
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From  this,  it  is  clear  that  F(z)  and  F  are  not  functions  of  x  and  y  and 
can  be  brought  outside  the  "del"  operator  allowing 


\V'  =  ^9-    (F(z)  -  F)  Ik  x  V  T'  (x,y,0)  .  (30) 

f  H 


Consequently, 


\V'»V  T'  =  ^-(F(z)-F)  [  Ik  x  V  T'  (x,y,0)-V  T'  (x,y,0)  (1  +  J)  ] 

=  ^-(F(z)-F)  (1+  %)  [   Ik  x  V  T'  (x,y,0)'V  T'  (x,y,0)]  (31) 


and  since  term  (a)  of  (31)  is  parallel  to  the  anomalous  temperature 
field  and  term  (b)  of  (31)  is  perpendicular  to  the  same  field,  then 

*  /v 

W»»V  T"  =  0  c   Since  \y'  =0  ,  we  have 
H 

\V  .  VHT'  =  0  , 

which  allows  u'-r—   and  v'-k—   to  be  neglected  in  (26).   Figures  14 
dx         dy 

and  15  clearly  show  that   initially  the  anomalous  currents  are  geostrophic 
and  parallel  to  the  anomalous  isotherms.   It  should  be  noted  that  these 
currents  need  not  remain  parallel  to  the  anomalous  isotherms  or  geos- 
trophic during  the  prediction  period;  however,  careful  examination  of 
the  predictions  indicates  they  do.   It  is  not  surprising  that  this  type 
of  motion  remains  geostrophic,  but  it  is  surprising  that  it  remains  equi- 
valent barotrophic.   Therefore,  for  the  purpose  of  this  interpretation, 

we  will  assume  that  W *  "  V  T'   is  negligible  during  the  forecast  period. 

H 
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3t  ' 
Additionally,   w'-k—   can  be  further  evaluated.   From  the  continuity 

oz 

equation  (8) , 


H  +  v„  •  *  "  °  • 


then 


t1  +  v„  *  w  +  It  +  v»  *  W '  =  o 

dz     H        dz      H 


where 


implying 

3w' 

If  +  VH  •  W  -  °  • 

Thus  integrating  up  from  the  level  z  to  the  surface  z  =  0  , 

o 
w'  =  J       V   •  w*  dz  +  w1 (0)  (32) 

z 

where  w* (0)  =  0  as  defined  by  the  boundary  conditions.   Therefore,  (32) 
can  be  written 


V  =  j   Vh  •  W'  dz  ,  (33) 

z 

/s  * 

and  since  V   •  W  =  0  ,  V   •  \u'  =  V   •  \u '   allowing  the  substitution  of 
H  H  H 

(30)  into  (33)  as 

o 
w'  =  ag  (  V   •  [(F(^~F)  ik  x  V  T'(x,y,0)]  dz  .  (34) 
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Since  the  coriolis  parameter  (f)  is  a  function  of  y,  and  T' (x,y,0)  is 
not  a  function  of  z,  Eq.  (34)  becomes 

V  =  Vu  •  £2- Ik  x  V„  T'(x,y,0)  F(z)]  (35) 

where 

o 
F(z)  =  /   (F(z)  -  F)  dz  . 


/ 


z 

Hence,  evaluating  (35) 


,■  -  [7  ,.a3TMx,y.O)     ^3TMx,y,0)     - 
H   f    ay  f     ox     w 


If  we  denote  the  geos trophic  current,  integrated  from  the  level  z  to  the 
surface  z  =  0  by  (Ug,  Vg)  , 


and 


vg(2)  .   at  ar- (x,y,0)  -M  _ 


then  (36)  can  be  written 


W  =  [-^  (Ug)  +  —  (Vg)]  .  (37) 


Therefore, 


and  so, 


-■J|afrattJ«-f'»« 


w'  3-  =  -  -  Vg(z)  -r—   .  (38) 

dz       f         dz 


27 


The  approximation  for  the  vertical  velocity  w'  applies  after  integration 
begins  as  well  as  initially.   This  result  shows  that  w'  is  the  vertical 
motion  produced  by  the  divergence  of  the  goestrophic  current.   Thus,  a 
northward  moving  layer  of  surface  water  (Vg  >  0)  converges  to  produce 
sinking  motion  below  the  surface  (w1  <  0) ,  while  southward  moving  surface 
water  (Vg  <  0)  diverges  and  produces  rising  motion  below  the  surface 
(w1  >  0) .   As  a  result  of  the  above  evaluated  terms  of  (26),  the  final 
prediction  equation  for  the  anomalous  temperatures  becomes 
12       3       4      5  6 

3t'    r-  3t*  -  3t*   ,  3t    .  3T   -  3t'   0  „  ,  ,  3t 

-~-  =  -[u  T-+v^—  +u  T~   +v  T~  +w-r—  -  —  Vg(z)  ^~ 
3t        3x     3y      3x      dy     3z    f       dz 

f        dz        ~  2     .2  2 

dx      dy         3z 


The  first  two  terms  on  the  right-hand  side  of  Eq.  (39)  represent 
the  horizontal  advection  of  the  anomalous  temperature  by  the  normal 
current.   These  are  the  only  terms  used  by  Namias  (1972a)  to  make  his 
prediction.   Terms  3  and  4  are  the  horizontal  advection  of  the  normal 
temperature  by  the  anomalous  current.   Term  5  is  the  vertical  advection 
of  the  anomalous  temperature  by  the  normal  current.   The  vertical  advec- 
tion of  the  normal  temperature  by  the  anomalous  current  and  the  vertical 
advection  of  the  anomalous  temperature  by  the  anomalous  current  are 
terms  6  and  7  respectively.   Term  8  is  the  horizontal  eddy  diffusion  of 
the  anomalous  temperature  and  term  9  is  the  vertical  eddy  diffusion  of 
the  anomalous  temperature  which  also  includes  the  anomalous  part  of  the 
surface  heat  flux  (see  discussion  below  (15)). 
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The  motivation  behind  the  inclusion  of  all  the  terms  of  (39)  is,  of 
course,  to  make  a  more  accurate  prediction,,   As  a  result  of  Jacob's  (1967) 
work,  which  showed  an  improvement  over  Namias1  predictions  (197  2a)  by 
including  terms  3  and  4  of  (39) ,  and  since  vertical  advection  is  so  im- 
portant in  the  theory  of  the  permanent  thermocline,  the  implication  is 
that  all  the  advection  terms  should  be  included  when  making  a  prediction 
of  the  sea-surface  temperature  anomaly. 
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III.   RESULTS 


Before  discussing  the  results  of  the  predictions  generated  by  the 
complete  dynamical  model,  a  brief  description  of  the  observed  anomalous 
temperature  fields  and  the  results  of  Namias ' (1972a)  predictions  fol- 
lows.  The  initial,  30,  60,  and  90  day  observed  anomalous  temperature 
fields,  as  seen  in  Figs.  15,  16,  17  and  18,  are  averages  over  three 
months  centered  about  the  months  of  October,  November,  December,  and 
January  respectively.   That  is,  the  initial  anomalous  temperature  field 
is  the  average  of  the  anomalous  temperatures  for  the  months  September, 
October,  and  November  (i.e.,  the  Fall  season)  of  1971. 

From  Figs.  15  through  18  two  cold  anomalies  and  one  warm  anomaly 
appear  to  dominate  the  area  of  interest.   The  cold  anomaly  near  the 
western  boundary  is  observed  to  decrease  in  amplitude  almost  linearly 
during  the  first  60  days  with  a  small  southward  movement.   The  double 
cell  cold  anomaly  near  the  northeast  boundary  decreases  in  amplitude 
in  the  western  part  but  maintains  its  relative  amplitude  and  position  in 
the  eastern  part.   Finally,  the  warm  anomaly  near  the  center  of  Fig.  15, 
increases  in  amplitude  with  the  leading  edge  of  the  maximum  anomalous 
temperature  contour  moving  eastward  some  13°  longitude. 

If  one  observes  Namias'  90  day  prediction  of  the  warm  anomaly  (Fig. 
4),  as  described  above,  the  easterly  movement  of  the  anomaly,  as 
advected  by  the  normal  winter  sea-surface  currents  (Fig.  3),  is  a  good 
approximation.   However,  the  amplitude  of  the  anomaly  in  his  prediction 
remains  relatively  constant,  when  it  was  observed  to  increase.   In  an 
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attempt  to  explain  this  discrepancy,  Namias  calculated  the  surface  heat 
fluxes,  but  the  results  obtained  indicated  a  negligible  (and  at  times 
opposed)  role  of  these  fluxes.   It  is  therefore  suggested  that  the  over- 
simplified use  of  the  prediction  equation,  as  described  earlier,  and  the 
lack  of  data  at  depth  is  responsible  for  not  predicting  accurately  the 
amplitude  of  the  anomaly. 

Attempting  to  improve  the  predictions  made  by  Namias  (1972a) ,  the 
completely  dynamical  model  of  this  thesis  was  run  for  a  corresponding 
90-day  period  using  the  initial  anomalous  temperature  field  (Fig.  15) 
and  the  model  generated  climatological  data  along  with  the  boundary  con- 
ditions as  described  earlier.   Four  cases  using  variations  in  the 
boundary  conditions  and  this  prediction  equation  are  discussed  below. 

A  brief  description  of  the  four  cases  will  be  given  first,  followed 
by  a  qualitative  discussion  of  the  results  obtained. 

CASE  1: 

In  CASE  1  the  anomalous  temperature  prediction  contained  all  the 
terms  implied  by  Eq.  (39).   The  results  of  CASE  1  are  shown  in  Figs.  19- 
24. 

CASE  2: 

In  CASE  2  the  model  calculated  the  90  day  prediction  as  in  CASE  1 
except  in  the  surface  thermal  boundary  condition,  Eqs.  (12),  (13),  and 
(14) ,  the  normal  climatological  sea-surface  temperatures  were  substituted 
for  the  predicted  sea-surface  temperatures.   This  effectively  removed  the 
anomalous  part  of  the  surface  heat  flux.   The  results  of  CASE  2  are  shown 
in  Figs.  25-30. 
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CASE  3: 

CASE  3  is  similar  to  CASE  2  except  normal  climato logical  temperatures 
were  substituted  in  place  of  the  predicted  temperatures  in  the  diffusion 

terms  of  (9) .   This  effectively  eliminates  the  anomalous  diffusion  terms, 

2      2        2 
d  T'     d  T'       d  T ' 
[AO2-^-  +  ■2— r-  ),  K  — j  ]  f  from  (39).   The  results  of  CASE  3  are  shown  in 

3x     8y       dz 

Figs.  31-36. 

CASE  4: 

CASE  4  is  similar  to  CASE  3  except  normal  climatological  currents 

(u,v)  were  substituted  for  the  predicted  currents  in  the  diffusion  terms 
2  2 

(A  V  u  ,  K  j  and  A  V  v  ,    <   j)  of  the  primitive  Eqs.  (5)  and  (6). 

dz  dz 

This  has  the  effect  of  removing  the  anomalous  diffusion  terms  (A  V  u'/ 
2  2 

K  — -  and  A  V  v',  K  r  )  from  the  equations  governing  the  anomalous 

9z  dz 

currents.   The  results  of  CASE  4  are  shown  in  Figs.  37-42. 

Concentrating  on  the  results  of  the  anomalous  temperature  predictions 
of  CASE  1  (Figs.  15,  20,  22,  and  24) ,  one  observes  that  all  the  anomalies 
appear  to  decrease  in  amplitude  with  the  cold  anomalies  disappearing  com- 
pletely by  the  end  of  that  90  day  prediction.   The  amplitude  of  the  warm 
anomaly  initially  was  found  to  be  .6°C  located  at  168°W,  40°N  and  is 
reduced  to  .1°C  after  90  days  with  its  center  at  162°W,  37°N.   As  a 
result  of  the  above  observations,  cases  2,  3,  and  4  were  designed  to 
determine  the  cause  of  the  amplitude  decrease  or  disappearance  of  these 
anomalies.   Furthermore,  the  western  edge  of  the  warm  anomaly  is  observed 
to  move  eastward  while  the  eastern  edge  remains  relatively  stationary. 
Below  20°N  some  anomalous  activity  appears  to  have  been  generated; 
however,  since  this  lies  outside  the  area  of  initial  observed  data,  it 
is  disregarded  and  further  discussion  of  this  region  will  be  omitted. 
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As  in  CASE  1/  the  cold  anomaly  at  the  northeast  boundary  disappears 
in  CASE  2  (Figs.  15,  26,  28,  and  30) .   Analyzing  the  initial  vertical 
profile  of  the  total  temperature  field,  it  was  found  that  the  lapse  rate 
was  unstable  thereby  implying  that  convective  adjustment  (term  6  (T)  in 
(9)),  is  the  cause  of  the  cold  anomaly  vanishing.   This  rapid  disappear- 
ance of  the  cold  anomaly  immediately  raises  the  question  of  the  impor- 
tance of  salinity.   If  salinity  were  included  in  the  present  model,  it 
is  speculated  that  the  cold  anomaly  could  maintain  itself.   That  is,  the 
cold  anomaly  at  the  surface  would  not  necessarily  be  unstable  if  a  (nega- 
tive) salinity  anomaly  were  present  to  stabilize  the  water  in  the  column. 
Another  possibility  is  that  salinity  is  negligible,  but  the  cold  anomaly 
extends  all  the  way  to  the  ocean  floor.   However,  this  explanation  seems 
less  likely  since  it  indicates  the  presence  of  such  a  large  anomaly  of 
heat  content. 

Investigation  of  the  initial  vertical  temperature  profile  of  the 
total  field  in  the  vicinity  of  the  cold  anomaly  near  the  western  boundary 
reveals  that  it  too  is  convectively  unstable.   However,  this  anomaly  damps 
less  rapidly  in  CASE  2  than  does  the  cold  anomaly  at  the  northeast 
boundary  which  would  indicate  its  decrease  in  amplitude  is  not  entirely 
due  to  convective  adjustment.   Additionally,  the  amplitude  of  the  central 
warm  anomaly  was  found  to  decrease  less  (to  a  value  of  ,2°C)  than  in 
CASE  1.   Except  for  the  above  two  differences,  the  60  and  90  day  predic- 
tions of  cases  1  and  2  are  relatively  similar.   Therefore,  with  respect 
to  the  observations  of  the  first  two  cases,  it  is  concluded  that  the 
anomalous  surface  heat  flux  induced  by  the  anomaly  itself  is  not  the 
major  cause  of  the  damping  of  the  anomalies,  and  that  convective  adjust- 
ment is  the  major  cause  of  the  cold  anomalies  vanishing.   The  effects  of 
diffusion  will  be  discussed  in  cases  3  and  4  below. 
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Figures  15,  32,  34,  and  36  show  the  results  of  CASE  3,  in  which  there 
is  no  eddy  diffusion  of  the  temperature  anomaly.   It  is  seen  that  the 
cold  anomaly  at  the  northeast  boundary  still  disappears  as  in  cases  1  and 
2.   The  western  cold  anomaly,  however,  appears  to  maintain  its  magnitude 
for  a  greater  duration  than  in  the  previous  two  cases;  the  greatest  dif- 
ference occurring  during  the  60  to  90  day  prediction  period.   The  western 
edge  of  the  warm  anomaly  moves  eastward  about  the  same  distance  as  before 
but  its  shape  is  altered  and  the  eastern  edge  remains  unchanged  as  in 
cases  1  and  2.   The  major  difference  observed  in  CASE  3  is  that  the  maxi- 
mum amplitude  center  of  the  warm  anomaly  only  decreases  from  ,6°C  ini- 
tially to  ,4°C  at  the  end  of  the  90  day  prediction.   Consequently,  the 
results  of  CASE  3  imply  that  the  anomalous  diffusion  terms  of  (39)  were 
significant  in  producing  the  heavy  damping  in  cases  1  and  2.   Since  re- 
moving the  anomalous  diffusion  terms  of  (39)  improved  the  results  of  the 
prediction  considerably,  but  did  not  completely  eliminate  all  of  the 
damping,  CASE  4  was  devised  in  which  the  diffusion  of  the  anomalous  cur- 
rents was  also  removed.   The  removal  of  these  additional  diffusion  terms 
should  eliminate  some  damping  because  of  the  dependence  of  the  tempera- 
ture on  the  currents  through  the  geostrophic  relationship  as  discussed 
previously. 

With  respect  to  dissipation,  movement,  and  amplitude  of  the  anomalies, 
CASE  4  (Figs.  15,  38,  40,  and  42)  produced  results  nearly  identical  to 
those  of  CASE  3.   Therefore,  it  is  concluded  that  the  diffusion  terms 
governing  anomalous  currents  in  the  primitive  equations  are  seemingly 
unimportant  in  the  prediction  of  sea-surface  temperature  anomalies.   Never- 
theless, CASE  4  represents  a  more  consistent  model  in  that  there  is  no 
diffusion  of  either  the  anomalous  currents  or  temperature. 
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If  one  now  compares  the  results  of  Namias'  prediction  (Fig.  4)  with 
the  results  found  in  this  study  using  a  completely  dynamical  model,  it 
is  first  observed  that  the  0°  isotherm  of  the  warm  anomaly  has  about  the 
same  configuration  and  location.   Secondly,  the  location  of  the  center 
of  maximum  amplitude  in  the  predicted  warm  anomaly  is  much  closer  to  the 
observed  90  day  position  in  Namias1  study  than  in  this  study.   However, 
if  one  closely  examines  the  location  of  this  anomaly  center  in  Figs.  15 
and  16,  it  is  seen  that  the  anomaly  center  moves  eastward  almost  15°  of 
longitude  during  the  first  30  days  and  then  remains  relatively  station- 
ary.  This  observation  would  imply  that  the  movement  of  the  anomaly's 
center  is  not  due  simply  to  advection  alone  as  indicated  in  Namias1  pre- 
diction.  Lastly,  the  amplitude  of  the  predicted  warm  anomaly  decreases 
in  all  cases  of  this  study  but  the  amplitude  of  the  same  anomaly  in 
Namias1  prediction  remains  constant  as  expected  when  employing  horizontal 
advection  alone.   Therefore,  the  vertical  advection  in  this  model  plays 
an  important  role  in  the  evolution  of  such  an  anomaly.   Whether  that 
role  in  the  present  prediction  or  any  future  predictions  is  simulated 
accurately  by  the  model  depends  critically  on  the  vertical  profile 
assumed  for  the  anomaly. 
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IV.   CONCLUSIONS 

This  study  has  shown  that  the  method  of  prediction  using  model  gener- 
ated climatology  along  with  observed  temperature  anomalies  in  a  fully 
dynamical  numerical  model  of  the  ocean  circulation  can  be  successful. 
Also,  the  simple  geostrophic  initial  conditions  for  the  anomalous  cur- 
rents is  successful  away  from  the  equator.   However,  the  quality  of  the 
above  results  depends  on  the  precise  initial  data  available. 

In  view  of  this  study,  the  need  for  temperature  data  (normal  and 
anomalous) in  the  mixed  layers  of  the  ocean  is  concluded  to  be  the  dominat- 
ing factor  influencing  the  prediction  of  the  sea-surface  temperature 
anomalies.   This  appears  quite  evident  from  our  observation  discussed 
above  in  reference  to  Figs.  15  and  16.   Since  this  rapid  movement  of 
the  anomaly  center  cannot  be  explained  by  horizontal  advection  alone, 
it  is  speculated  that  a  sub-surface  warm  anomaly  may  have  existed  down- 
stream which  was  brought  to  the  surface.   Thus,  vertical  advection  may 
have  played  a  significant  role;  however,  the  lack  of  data  beneath  the 
surface  makes  these  speculations  virtually  impossible  to  investigate. 

The  present  model,  with  no  salinity,  cannot  accurately  predict  the 
evolution  of  high  latitude  cold  anomalies  which  do  not  extend  to  the 
ocean  floor.   Therefore,  initialization  should  include  salinity;  however, 
at  present,  large  scale  anomalous  salinity  data  is  not  available. 

When  comparing  the  results  of  cases  3  and  4,  it  is  observed  that 
the  removal  of  the  diffusion  terms  governing  the  anomalous  currents  had 
little  effect,  if  any,  on  the  90  day  prediction.   Thus,  further 


36 


predictions  using  a  dynamical  model  should,  for  consistency,  omit  all  of 
the  anomalous  diffusion  terms  (CASE  4) . 

Although  the  results  of  this  study  did  not  produce  predictions  which 
corresponded  exactly  with  the  observed  anomaly  field,  the  use  of  a  fully 
dynamical  model  in  making  these  predictions  has  given  us  a  greater  in- 
sight as  to  the  role  played  by  vertical  advection  and  diffusion.   It  is 
therefore  believed  that  further  studies  should  be  conducted  using  this 
fully  dynamical  model  with  more  accurate  initial  data  and  with  salinity 
included. 

Additionally,  this  study  uses  only  climatological  data  to  define  the 
atmospheric  forcing.   However,  an  anomalous  atmosphere  would  produce 
anomalous  heating  and  wind  stresses  which  could  produce  anomalous  features 
in  the  ocean.   For  example,  Namias  (1965)  investigated  the  effects  of  an 
anomalous  wind  stress  and  showed  that  this  anomalous  atmospheric  forcing 
may  be  important  in  some  cases.   Thus,  further  investigations  should  be 
conducted  to  determine  the  response  of  the  ocean  to  various  anomalous 
atmospheric  forcing. 
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